Pseudogap enhancement due to magnetic impurities in d-density waves 
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We study the effect of quantum magnetic impurities on d-wave spin density waves (d-SDW). 
The impurity spins are aligned coherently according to the spin space anisotropy of the condensate. 
Both the order parameter and transition temperature increases due to the coherent interplay between 
magnetic scatterers and d-SDW. This can explain the recent experimental data on the pseudogap 
enhancement of Ni substituted NdBa2{Cui_yNij,}06.8 from Pimenov et al. (Phys. Rev. Lett. 94, 
227003 (2005)). 
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The nature of the pseudogap phase of high Tc super- 
conductors is highly controversial. Among many oth- 
ers, Chakravarty et ali proposed that it results from d- 
density wave order. This phase is a charge density wave 
(CDW)^ with a gap of d-wave symmetry. This causes 
the lack of periodic charge modulation, and evokes the 
notion "hidden-order". However, there are many other 
unconventional density wave phases^, with similar prop- 
erties to d-density waves^*^. For example, d-wave spin 
density waves (d-SDW) are almost indentical to d-CDW 
except the hidden spin space anisotropy, and can account 
for several experimental results in the pseudogap phase 
of high Tc superconductors equally well as d-CDW does. 

Last year, Pimenov et alii reported that the pseudogap 
energy in underdopcd NdBa2Cu06.8 is enhanced by sub- 
stituting Cu in the Cu02 plane with Ni by measuring the 
c-axis optical conductivity. We have established recently 
that the pseudogap phase should be d-wave density wave 
through the analysis of both the giant Nernst effect and 
the angle dependent magnetoresistanca^. However, we 
cannot decide whether it is d-wave charge or spin den- 
sity wave. The earlier Pauli limiting behaviour of d- 
density waves suggested d-CDWii2*ii*i^. However, it is 
easily seen that d-SDW also exhibits Pauli limiting when 
the spin anisotropy axis in d-SDW lies in the a — 6 plane. 
Further a recent analysis of impurity scatteringi^ under 
general conditions suggests that such an enhancement of 
d-density wave order is impossible, if the condensate is 
d-CDW and/or the impurity is nonmagnetioi^*i^. This 
leaves us the unique possibility that d-density wave needs 
to be d-SDW and the impurity scattering is of Kondo 
type. 

Therefore we consider in the foUowings d-SDW in the 
presence of magnetic impurities with the Kondo coupling. 
The study of magnetic imputities has a long historjiiSiii. 
In order to couple the impurity spins to the underlying 
d-SDW, the Kondo coupling has to be non-local and has 
to have a d-wave component as in Ref. m The interac- 
tion between quantum magnetic impurities and d-SDW 
is similar in essence to the RKKY interaction, but the 
electrons mediating the coupling between the spins are 
not free but participate in collective phenomenon, and 



belong to the density wave condensate. This leads to 
the enhancement of the pseudogap energy, similarly to 
ferromagnetic superconductorsiS. There, the Pauli term 
due to an external magnetic field can compensate the ex- 
change term due to ferromagnetism, and superconductiv- 
ity will be revived. It is worth mentioning that d-density 
waves have unusual magnetic properties with the inclu- 
sion of spin-orbit coupling as well2S. 

In general, nonmagnetic impurities such as Zn, lead 
to the destruction of the underlying phaseSiiS^, regard- 
less to its symmetry, and are thought to belong to the 
unitary scattering limit. Ni impurities have also simi- 
lar impact on superconductors due to the lack of spin 
space anisotropy, and act as simple potential scatterers 
in the Born limit. In the situation considered below, 
their magnetism can couple directly to the order param- 
eter of d-SDW, completely modifying the simple scatter- 
ing picture. No such mechanism is expected to exist in 
superconductors, hence magnetic Ni impurities suppress 
superconductivity similarly to Zn. 

As a model, wc consider an effective low-energy 
Hamiltonian describing d-wave spin density waves (d- 
SDW)2i24, interacting with randomly distributed quan- 
tum magnetic impurities: 
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Here e(k) is the kinetic energy spectrum, A is the d- 
SDW order parameter, /(k) — {cos,{kxa) — cos(fcj,a))/2 
is the gap function with d-wave symmetry. For simplic- 
ity, and to make connection with earlier analysisi, we 
study a commensurate situation with best nesting vec- 
tor Q = (tt/a, 7r/a). A magnetic impurity (Sj) at Rj 
interacts with the electron spin density (s(r)) through 
the non-local Heisenberg exchange coupling. Its matrix 
elementi^ describing scattering with Q is expanded in 
terms of Fermi surface harmonics, and we retain the term 
possessing d-wave symmetry, namely J(k, Q) oc /(k). 



2 



the others are unable to couple to the wavevector depen- 
dence of the gap. To hivestigate the effect of magnetic 
interaction on d-SDW, we study its change on the grand 
canonical potential perturbatively. To lowest order, one 
finds 
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which suggests that each impurity feels a local aligning 
"magnetic" field manifested through A/P cos(QRj), 
is the component of the impurity spin parallel to the 
spin of underlying d-SDW which lies in the a — b plane. 
Depending on the sign of JA, the impurity will order lo- 
cally in ferro- or antiferromagnetic fashion with respect 
to the local fieldi^, but the overall arrangement of impu- 
rities will follow an antiferromagnetic pattern due to the 
cos(QRj) factor, as is visualized schematically in Fig. ^ 




FIG. 1: (Color online) The impurity spin configuration is 
sketched in the ordered phase. The solid line denotes the 
alternating field dictated by the impurity-density wave inter- 
action, the colored objects represent the impurity spins. 

To make this picture more quantitative, we perform 
a two order parameter mean field theory, one for the d- 
SDW and another for the impurities by assuming, that 
the impurity spins are aligned as = tocos(QRj ). The 
mean field decoupled magnetic interaction reads as 
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Here is the number of unit cells, P is the d-SDW 
interaction^ and Ni is the number of impurity atoms. To 
average over impurities, we use the standard approach to 
consider only non-crossing ladder type diagrams^^. Usu- 
ally, in the first order Born approximation, nonmagnetic 
impurities simply shift the chemical potentialM. In our 
case, however, magnetic impurities couple to the order 
parameter. A, and in the first order Born approximation, 
the new energy spectrum is given by 
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The use of the Born approximation is justified from the 
fact that Ni impurities in unconventional superconduc- 



tors can convincingly be described with it, while Zn sub- 
stituted samples call for the unitary limit2iiS^. The sta- 
bility conditions of the order parameters can be obtained 
from the grand canonical potential, which takes the from 

k e RBZ 

Here S is the impurity spin quantum number, RBZ 
stands for the reduced (antiferromagnetic) Brillouin zone. 
The first two terms are characteristic to density waves^, 
the third one stands for the electron-impurity interaction, 
while the latter describes to certain extent the quantum 
nature of the magnetic impurity, as opposed to classical 
spin9ii. After minimization with respect to A and m, 
the coupled gap equations are obtained as 
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TT-i = Ni/N, f{x) is the Fermi distribution function. The 
first one is the usual BCS type gap equation, while the 
second one with Brillouin function on the right hand side 
accounts for the quantum nature of the magnetic impu- 
rity. Regardless to the detailed form of the kinetic energy 
spectrum, Eq. Q can be solved at T = to give m = S, 
the impurity spin is aligned completely. The above gap 
equations can further be simplified in the continuum 
limi1i& and assuming perfect nesting (e(k) -f e(k — Q) = 
0) , but we would like to emphasize that our result are ro- 
bust with respect to variation of the kinetic energy spec- 
trum, commensurability and imperfect nesting. With 
this, the zero temperature d-SDW order parameter is de- 
termined from 
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with pq the normal state density of states, Aq is the or- 
der parameter of the pure system at T = 0. In spite of 
its relative simplicity, this formula constitutes the main 
result of this paper and is shown in Fig. |21 for various 
interactions. For any finite J, A is enhanced with re- 
spect to its pure value Aq. The steady increase stems 
from Eq. Q, since impurities act like a source term. 
These findings are in accordance with Ref. 0, where a 
steady increase of the pseudogap has been reported with 
respect to Ni impurity concentration, by focusing on the 
suppression of the electric conductivity. Our result indi- 
cates that this can arise from the coherent interplay of 
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magnetic impurities and d-SDW. Moreover, the experi- 
mentally measured pseudogap energy scale is not simply 
A, but A + iiiJS as seen from Eq. (0J, causing further 
enhancement of the pseudogap. As opposed to this, Zn 
impurities act like potential scatterers, and provide us 
with pair-breaking effect, leading to the destruction of 
the superconducting or density wave condensate. 
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2: (Color online) The concentration and magnetic cou- 



pling dependence of the order parameters is shown (A with 
blue solid, m with red dashed line) at T = in the weak- 
coupling limit for Ppo — 0.7, 0.5 and 0.3 from bottom to top. 
Note the strong increase in A even for low impurity concen- 
trations. 

The change in the transition temperature (Tc) is dif- 
ferent from the Abrikosov-Gor'kov formula deduced for 
non-magnetic impurities, and reads as 

1 2n,J^S{S + l) 



2n,J^S{S + l)' 

The right hand side of this equation is always posi- 
tive, hence similarly to A, Tc also enhances in the pres- 
ence of magnetic impurities. For arbitrary temperatures, 
the gap-equations have been solved numerically, and are 
shown in Fig. for S ~ 1, characteristic to Ni. When 
the Hcisenberg coupling (J) is large, the spins retain their 
maximum value (S) up until the close vicinity of Tc, and 
a small amount of impurities increase significantly A. On 
the other hand, for weaker J, the d-SDW order parame- 
ter follows the usual BCS like temperature dependence, 
and m reaches its maximum value only at low tempera- 
tures. 

We mention here, that the finite value of m results in 
small amplitude antifcrromagnetic modulation of the lo- 
calized spins. Recently such weak low temperature stag- 
gered magnetization lying in the a — b plane has been 
observed in YBa2Cu306.3^ with an amplitude toq = 
0.05/iB well above the superconducting transition tem- 
perature (55 K) at 310 K. We believe that this antifcr- 
romagnetism is related to the interaction between the 
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3: (Color online) The temperature dependence of the 



order parameters is shown (A with blue solid, m with red 
dashed line) in the weak-coupling limit for Ppo = 0.5 and 
S — 1. The top panel shows weak impurity-electron coupling 
Ao/J — 4.5 with m = 0, 0.03, 0.1 and 0.2 with increasing 
Tc- In the bottom panel, strong impurity-electron coupling is 
used with Ao/J = 0.45, rn = 0, 0.01, 0.03 and 0.06 with in- 
creasing transition temperature. For stronger coupling, small 
concentration causes dramatic increase in both A and Tc- 



pseudogap phase and localized moments. Indeed, the 
measured magnetic intensity closely resembles to that in 
Fig. 121 Albeit orbital antiferromagnetismi does not exist 
in our d-SDW model as opposed to its d-CDW counter- 
part, the sensitivity of its order parameter to local mag- 
netic fields provides us with weak antiferromagnctism, 
similarly to URu2Si2i^. 

Similar phenomenon is expected to occur in itiner- 
ant antiferromagnets (namely conventional spin density 
waves with constant order parameter—) interacting with 
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magnetic impurities. The impurities will align follow- 
ing the spin density oscillations of electrons, leading to 
an enhancement of the ordering amplitude. We pre- 
dict the realization of the aforementioned coherent in- 
terplay between magnetic impurities and condensate in 
the SDW phase of Bechgaard salts (TMTSF)2X with 
X=PF6, CIO4, AsFg etc. polluted with magnetic scat- 
terers. 

In conclusion, we have demonstrated that magnetic 
impurities are coherently aligned in a d-wave spin den- 
sity wave, causing the increase of the order param- 
eter of the condensate. This can be identified with 
the pseudogap energy scale, whose enhancement is in 
agreement with a recent experiment in Ni substituted 
NdBa2{Cui_yNiy}06.3i. As was pointed out, magnetic 
correlations play an important role in the explanation. 
We predict, that beyond the pseudogap enhancement in 



NdBa2{Cui_j,Nij,}06.8: small amplitude antiferromag- 
netic ordering should be observable due to the Ni im- 
purities, similarly to YBa2Cu306.3^. 
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